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ABSTRACT 


Cosmological Gamma-Ray Bursts (GRBs) are known to arise from distinct progenitor channels: short GRBs mostly 
from neutron star mergers and long GRBs from a rare type of core-collapse supernova (CCSN) called collapsars. 
Highly magnetized neutron stars called magnetars also generate energetic, short-duration gamma-ray transients called 
Magnetar Giant Flares (MGFs). Three have been observed from the Milky Way and its satellite galaxies and they 
have long been suspected to contribute a third class of extragalactic GRBs. We report the unambiguous identification 
of a distinct population of 4 local («5 Mpc) short GRBs, adding GRB 070222 to previously discussed events. While 
identified solely based on alignment to nearby star-forming galaxies, their rise time and isotropic energy release are 
independently inconsistent with the larger short GRB population at 299.996 confidence. These properties, the host 
galaxies, and non-detection in gravitational waves all point to an extragalactic MGF origin. Despite the small sample, 
the inferred volumetric rates for events above 4 x 10*^ erg of Ruar = 3.8139 x 10° Gpc^? yr-! place MGFs as the 
dominant gamma-ray transient detected from extragalactic sources. As previously suggested, these rates imply that 
some magnetars produce multiple MGFs, providing a source of repeating GRBs. The rates and host galaxies favor 
common CCSN as key progenitors of magnetars. 


Keywords: gamma rays: general, methods: observation 


1. INTRODUCTION 


The history of GRBs and magnetars are intertwined. 
Short bursts of gamma-rays were recorded by the Vela 
satellites beginning in 1967 (Klebesadel et al. 1973), 
and were given the phenomenological name GRBs. 
GRB 790305B was localized by the InterPlanetary Net- 
work (IPN) to the Large Magellanic Cloud (Mazets 
et al. 1979; Evans et al. 1980). It was unique in being 
the brightest event seen at Earth, the prompt emission 
had a long-lasting, exponentially-decaying, periodic tail 
(Barat et al. 1979) and additional, weaker bursts were 
localized to the same source (Mazets et al. 1979). Im- 
mediately there were papers investigating if the main 
event shared a common origin with other GRBs (Mazets 
et al. 1982; Cline et al. 1980). It is now known to be the 
first signal identified from a magnetar. 

Key results on the nature of GRBs in the subsequent 
decades were often proven by population-level statistical 
analysis before direct “smoking-gun” proof. Perhaps the 
greatest debate was whether these events had a galac- 
tic or an extragalactic origin, with the latter initially 
disfavored as it would require intrinsic energetics be- 
yond anything previously known. Proof came first indi- 
rectly via statistical studies on the spatial distribution 
of GRBs (Meegan et al. 1992) and then directly from 
redshift measurements (Metzger et al. 1997). 

Studies of the prompt GRB emission provided strong 
evidence in favor of two populations (Kouveliotou et al. 
1993), with short and long GRBs traditionally sepa- 
rated at 2s as measured by the Too parameter. Long 
GRBs were tied to broad-line type Ic core-collapse su- 
pernovae called collapsars (Galama et al. 1998). The 
Neil Gehrels Swift Observatory (Swift) mission enabled 
successful detections of afterglow from a sample of short 
GRBs. Circumstantial evidence pointed towards a neu- 
tron star merger origin (Eichler et al. 1989; Fong et al. 
2015) with direct confirmation that some GRBs arise 
from binary neutron star mergers came with GW170817 
and GRB 170817A (Abbott et al. 2017). 

Yet another debate on the behavior of GRBs is 
whether or not the sources repeated. This is best 
explained using modern parlance. Soft Gamma-ray 
Repeaters (SGRs) are galactic magnetars named phe- 
nomenologically for the weak, recurrent short bursts 
that first identified them before their physical origin 
was known. SGR flares are classified as distinct from 
GRBs, and have recently been tied to radio emission 
similar to the cosmological Fast Radio Bursts (Boch- 
enek et al. 2020). The flare on March 5, 1979 and the 
subsequent similar events GRB 980827 (Mazets et al. 
1999b; Hurley et al. 19993) and GRB 041227 (Palmer 
et al. 2005; Frederiks et al. 2007a) from magnetars in 


the Milky Way are referred to as Magnetar Giant Flares 
(MGFs). The designation for the prompt emission of 
MGFs often carries the GRB designation, which we use 
here. GRBs are now not thought to repeat as collapsars 
and neutron star mergers are cataclysmic events. While 
several galactic magnetars have been observed to pro- 
duce multiple SGR flares, none have been observed to 
produce multiple giant flares (though this is not surpris- 
ing). The historic debate on potential repeating GRBs 
was likely confounded by magnetar transients before the 
separation of SGR flares from GRBs. 

We here refer to GRBs 790305B, 980827, and 041227 
as the known MGF sample. The detection of three 
from the Milky Way and its satellite galaxies implies 
a high intrinsic rate on a per-galaxy or volumetric basis. 
'These events should be detectable to extragalactic dis- 
tances by GRB monitors such as Konus- Wind (Aptekar 
et al. 1995), Swift-BAT (Barthelmy et al. 2005), and 
Fermi-GBM (Meegan et al. 2009). However, at these 
distances only the immediate bright spike would be de- 
tectable and the event should resemble a short GRB 
(Hurley et al. 2005). There are two events discussed in 
previous literature as extragalatic MGF candidates, be- 
ing GRB 051103 (Ofek et al. 2006; Frederiks et al. 2007b; 
Hurley et al. 2010) and GRB 070201 (Mazets et al. 2008; 
Ofek et al. 2008), whose chance alignment coincidence 
was measured to be ~1% (Svinkin et al. 2015). 

''here have been population-level searches for addi- 
tional events, which identified no additional candidates 
(Popov & Stern 2006; Ofek 2007; Svinkin et al. 2015). 
However, these studies allow us to constrain the frac- 
tion of detected short GRBs that have an MGF origin: 
Ofek (2007) show that the rate of galactic events requires 
this to be >1%, while the lack of additional candidates 
found in several searches constrain the upper bound to 
be «896 (Tikhomirova et al. 2010; Svinkin et al. 2015; 
Mandhai et al. 2018). These studies and tehir conclu- 
sions generally assumed that the brightest MGFs could 
be detectable to tens of Mpc. 

Recently GRB 200415A was identified as the third 
and likeliest extragalactic MGF (Svinkin et al. 2021). 
In this work, we perform a new population-level search 
utilizing the largest GRB sample, new galaxy catalogs 
that are both more complete and provide additional in- 
formation, and develop a new formalism to determine if 
we can prove extragalactic MGFs contribute to the ob- 
served GRB population. Section2.4 details the search 
formalism which identifies four nearby events, identify- 
ing an additional extragalactic candidate. The progeni- 
tors of our identified sample are investigated in Section 3, 
the implications of which are discussed in Section 4. We 
conclude with discussions in Section 5. 
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2. LOCAL GRBS 


'The *smoking-gun" evidence of an MGF is the long 
periodic tails which are modulated by the rotation pe- 
riod of the neutron star (Hurley et al. 1999b) and also 
show quasi-periodic oscillations related to the modes of 
the neutron star itself (Barat et al. 1983; Strohmayer 
& Watts 2005; Israel et al. 2005; Watts & Strohmayer 
2006). However, these signatures are not unambiguously 
identifiable at extragalactic distances with existing in- 
struments. As such, we follow prior population-level 
searches and focus on spatial information: if a well- 
localized short GRB is an MGF it should occur within 
7-50 Mpc and be consistent with a cataloged galaxy. We 
combine existing GRB and galaxy catalogs to build the 
most complete set of information from existing litera- 
ture. For each individual burst we quantify our be- 
lief that it is an MGF from a known galaxy through 
comparison of two PDFs, which are discussed below. 
These PDFs are generated in HEALPix (Gorski et al. 
2005). The resolution of HEALPix maps is defined by 
the NSIDE parameter, where the number of total pixels 
is equal to the square of the NSIDE times twelve. The 
maps were generated with NSIDE=8192, corresponding 
to a pixel width of ~0.5 arcminutes. 


2.1. The GRB Sample 


We utilize data from CGRO-BATSE (Fishman et al. 
1989), Konus- Wind (Aptekar et al. 1995), Swift-BAT 
(Barthelmy et al. 2005), Fermi-GBM (Meegan et al. 
2009), and additional information from the IPN!. Trig- 
gers from the same events were matched utilizing tem- 
poral information for all events and spatial information 
(Ashton et al. 2018) when available. The total sample 
contains more than 11,000 GRBs observed, with 71,200 
short GRBs using the standard 2s cutoff. 

Our burst sample selection requires three things. We 
consider only short GRBs (Tyo < 2s) where the Too 
used is the shortest reported by any triggering in- 
strument. Second, we require the bolometric fluence 
(1keV-10 MeV) determined from a broadband instru- 
ment (Konus, BATSE, or GBM), converting from the 
instrument-specific ranges as necessary. Intercalibra- 
tion uncertainties are within 2596. For the trigger times, 
duration, and spectral properties we utilized the latest 
catalog information (Paciesas et al. 1999; Svinkin et al. 
2016; Lien et al. 2016; von Kienlin et al. 2020), updated 
online catalogs?, GCN circulars, and performed dedi- 
cated analysis when necessary. 


1 ss].berkeley.edu/ipn3/index.html 
2 http://www.ioffe.ru/LEA/shortGRBs/Current/index.html 


Lastly, we require well-localized GRBs, constructed 
from all available information. For BATSE localization 
we utilize the latest catalogs (Goldstein et al. 2013) and 
apply the largest systematic error (Briggs et al. 1999). 
Swift-B AT positions are taken from the updated Swift- 
BAT Catalog? and Swift-XRT localizations are utilized 
when available^. Fermi-GBM localizations are quasi- 
circular and were generated using the latest methods 
(Goldstein et al. 2020) for all bursts. 

KONUS localizations are an ecliptic band which are 
summarized in the IPN catalogs. The IPN compiles lo- 
calization information for GRBs, including the timing 
annuli derived from the relative arrival times of gamma- 
rays at distant spacecraft. Information used here is from 
the IPN localizations of Konus short GRBs through 2020 
(Pal'Shin et al. 2013) and the IPN list kept up to date 
online?. Additional IPN localizations were compiled for 
more than 100 additional short GRBs for this work, 
which were added to the online table. The location 
information, including systematic error, from the au- 
tonomous localizations, timing annuli, and Earth occul- 
tation selections are converted to the HEALpix format 
using the GBM Data Tools®. These independent PDFs 
are combined into a final PDF referred to as PCPP, 

The localization threshold is set to a 9096 confidence 
area « 4.125 deg? when including systematic error. This 
value is chosen as it is 1/10,000 the area of the sky, is 
comparable to the sum of the angular size of galaxies (as 
defined in the following section) within 200 Mpc, and is 
between previously used thresholds (Svinkin et al. 2015). 
With the bolometric fluence measure requirement and 
the removal of bursts with known redshift (Lien et al. 
2016) beyond the distance where the event may be a 
detected MGF, we are left with a sample of 250 short 
GRBs. We do not apply more stringent cuts on spectral 
or temporal information at this stage as the relevant 
parameters are not uniformly reported in GRB catalogs. 


2.2. The Galaxy Sample 


For the galaxies considered in this work we require the 
position (KA, Dec, Distance), angular extent (if non- 
negligible at our spatial resolution; represented here as 
ellipses), and the current Star Formation Rate (SFR). 
The z=0 Multiwavelength Galaxy Synthesis (z0MGS) 
Catalog (Leroy et al. 2019) combines the ultraviolet ob- 


3 
html 
4 https://swift.gsfc.nasa.gov/archive/grb table/ 


5 http://www.ssl.berkeley.edu/ipn3/ 
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servations from GALEX (Morrissey et al. 2007) with the 
infrared observations of WISE (Wright et al. 2010) to 
uniformly measure gas and dust for galaxies within ap- 
proximately 50 Mpc. As a result, for galaxies contained 
in this catalog these measures of the distance and SFR 
are our default values. The angular size of galaxies is 
represented as an ellipse when data allows or as a cir- 
cle when the axial ratio is not known. Angular extent 
is taken from the input catalogs, but is generally the 
Holmberg isophote, i.e. where the B band brightness is 
26.5 mag arcsecond?. 

The Census of the Local Universe (CLU) Catalog 
(Cook et al. 2019) aims to provide the most complete 
catalog of galaxies out to 200 Mpc. We use the CLU 
measures of distance and SFR when they are not pro- 
vided by z0MGS, and we use the CLU measures for 
angular size (which are not provided by the zOMGS). 
When missing, we add position angle information from 
HyperLEDA (Paturel et al. 2003). The SFR measures 
of these two catalogs correct for internal extinction us- 
ing WISE4/FUV luminosities. To ensure completeness 
within «10 Mpc we supplement these two catalogs with 
the Local Volume Galaxy (LVG) Catalog (Karachent- 
sev & Kaisina 2013). The three catalogs are matched 
by name, with help from the NASA/IPAC Extragalac- 
tic Database (NED)’, and position information. 

We consider galaxies between 0.5 Mpc (excluding the 
Milky Way and its satellite galaxies) and 200 Mpc (be- 
yond where MGFs can be detected), which leaves more 
than 100,000 galaxies. The SFR is a key parameter in 
our method and our inferences also rely on scaling the 
properties of our host galaxy. The Milky Way SFR used 
here is 1.6540.19 Mọo/yr (Licquia & Newman 2015). 
We specify the SFR for NGC 3256, which was identified 
in Popov & Stern (2006) as being a likely source of de- 
tectable extragalactic MGFs. We searched the literature 
for values of the active SFR in this galaxy and take the 
value of ^36 5/yr from Lehmer et al. (2015) which is 
inferred using UV information and is among the middle 
reported values. 


2.3. MGF Spatial Distribution 


We seek an all-sky PDF, PMGF, representing the 
probability that a given position is to produce a MGF 
with a particular fluence at Earth. Note that this is de- 
termined by the fluence of each burst considered, but is 
constructed independently of the location of the burst 
itself, PCAP, The comparison of the two PDFs gener- 
ated for each burst quantifies the likelihood that a given 
short GRB has an MGF origin, which is performed in 


7 https://ned.ipac.caltech.edu/ 


the next section. This section details the burst-specific 
construction of PMGF, 

If a given burst has an MGF origin it should arise from 
a cataloged galaxy and its intrinsic energetics should fall 
into the expected range. To construct this we compute 
a weight for each galaxy representing how likely it is to 
have produced the observed fluence for the burst under 
consideration. This weight has two-components: a lin- 
ear weighting with SFR and a more complex weighting 
that compares the inferred intrinsic energetics (deter- 
mined by the burst fluence and potential host galaxy 
distance) against an assumed PDF. 

Magnetars are expected to be able to produce MGFs 
only for a short period of time (approximately 10 kyr Be- 
niamini et al. 2019), tying the predicted rate of MGFs 
to the rate of their formation. The rate of CCSN can be 
inferred from the SFR since the lifetimes of stars that 
undergo core-collapse is much shorter than the timescale 
probed by the SFR tracers (Botticella et al. 2012). Un- 
der the assumption that the dominant formation channel 
for magnetars is CCSN (which is explored in Section 4) 
we can infer the rate of MGFs from a galaxy from its 
SFR. Thus, each galaxy is linearly weighted with SFR. 
We use the far ultraviolet measure of SFR (Lee et al. 
2010) when available as it should track massive stars 
likely to undergo core-collapse, otherwise we use the Ho 
measure (Kennicutt Jr 1998) scaled by the average dif- 
ference from galaxies with both measures to account for 
the lack of dust correction in the LVG catalog. 

Next we can determine the total isotropic-equivalent 
energetics of a potential burst-galaxy pair as Eis = 
And?S where S is the burst fluence and d the distance 
to the potential host. This value can be compared to 
an assumed intrinsic energetics PDF to determine how 
likely the event is to be an MGF. For example, a partic- 
ularly high fluence short GRB spatially aligned with a 
distant galaxy would require an intrinsic energetics far 
beyond what has been observed in the galactic MGFs, 
excluding an MGF origin. We note that some studies 
utilize the peak luminosity LM@* but we work with an 
Eiso distribution as there is stronger theoretical guid- 
ance on the maximum total energy that can be released 
(related to the magnetic fields of the magnetar) than on 
the timescale that it is released. 

We now construct an informed intrinsic energetics 
function, assuming a power-law distribution with an as- 
sumed minimum and maximum value, which is similar 
to the behavior of lower energy magnetar flares (Cheng 
et al. 1996). Our method bypasses the need for an as- 
sumed detection threshold, which is difficult to quantify 
when considering many instruments over 30 years. The 
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Figure 1. The initial assumed MGF energetics distribution, 
with Fisomin and Eiso max set to the x-axis boundaries. The 
PDF form is (1 — o) E? / (Eire — E452). As described in 
the text, a = 1.3 + 0.9 (at 90% confidence). The three Eiso 
values from the known MGFs used to constrain the slope are 
shown as black vertical lines. 


assumed and inferred values are reported below, with 
the initially determined distribution shown in Figure 1. 

The slope of a power-law can be determined via max- 
imum likelihood, independent of an assumed maximum 
value, as 


£ Figo | «0-1 | 
a(R) ve vou 


Eiso,min 
(1) 


where the sum is over the observed Eiso and Eiso min 
is the lowest considered value Newman (2005); Bauke 
(2007). We set Eisomin as 1.0 x 107? erg which is a 
factor of a few below the lowest value measured in a 
known MGF as shown in Table1 but above the bright- 
est SGR flare that lacked the periodic tail emission 
(Mazets et al. 1999a). Iterating over the Eiso values of 
the known MGFs (GRBs 790305B, 090827, and 041227) 
gives a = 1.3 + 0.9 at 90% confidence, where we have 
included the O(n~*) error contribution. In order to min- 
imize the required computation we assume the centroid 
(a = 1.3) in what follows; the effect of this assumption 
on our results is discussed in the closing paragraph of 
this section. 

There must be a physical maximum energy for an 
MGF, which should be related to the total magnetic 
energy. This is supported by the lack of detections 
of more energetic events otherwise consistent with an 
MGF origin. The highest EF; observed for a known 
MGF is 2.3 x 104° erg which comes from the magne- 
tar with the highest reported magnetic field at the 
surface of 2.0 x 1015G (Olausen & Kaspi 2014). We 
note this reported value is approximately 3 times larger 
than the dipolar spin-down inferred magnetic field value 


a=1l+n 


of 7x101^ G (Younes et al. 2017), but we have con- 
firmed this does not affect our results. To determine an 
Eiso,max for our search we assume a dipole field, where 
the available energy scales as B?, and a nominal maxi- 
mum magnetic field strength of —1.0x 1019 G. This gives 
Eis max = 2.3 x 10% erg x (L0 x 10/9G/2.0 x 1015G)? = 
5.75 x 10% erg. 

This allows us to determine the burst-specific two- 
component weight for each of the 100,000 galaxies in 
our sample, which are weighted linearly by its SFR mul- 
tiplied by the value of the Ei, PDF for the inferred 
energetics considering the burst fluence and galaxy dis- 
tance. The sum of the galaxy weights is normalized to 
unity. Then, PMGF is built by placing the calculated 
weights at the position of the host galaxy. If the angular 
diameter of the galaxy is larger than the effective reso- 
lution of our discrete sky representation (~arcminute?) 
then its weight is uniformly distributed over its angular 
extent. 


2.4. The Search 


For each of the 250 short GRBs in our sample we 
generate PRE from the observations of the GRB and 
PMGF from theoretically motivated expectations. We 
quantify the likelihood that a given GRB has an MGF 
origin using Q = 4m >), PP RP PMGF/A; where PCRB 
and PMC” indicate the probability for each PDF in the 
ith sky region, which has area A; (Ashton et al. 2018). 

Significance is determined by the empirical False 
Alarm method (e.g Messick et al. 2017) with Q as our 
ranking statistic. Our backgrounds are generated by 
simulating different galaxy distributions. Each iteration 
is generated by uniform rotation of the 2D (RA, Dec) 
positions of the galaxies in our sample, which main- 
tains the distance and SFR distributions as well as local 
structure. Population-level confidence intervals created 
through comparison of each rotation against our full 
GRB sample with results are shown in Figure2. At 
3 and 4 events the short GRB sample has an excess 
surpassing 50 discovery significance, with individual 
significance values of the four bursts between 1.2x1074 
and 4.9 x 107 as given in Table 1. 

Three of the four are discussed in the literature 
as extragalatic MGF candidates. The Konus- Wind 
lightcurves are shown in Figure 3. GRB 070201 has the 
least robust association to a nearby galaxy; however, 
the localization is comparatively large (~10x the other 
events) and M31 has the largest angular size of any 
galaxy in our sample, together lowering Q even for real 
associations. We confirm this by checking GRB 790305B 
with the Large Magellanic Cloud (Evans et al. 1980; 
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Figure 2. The discovery of a local but extragalactic popula- 
tion of GRBs. Q is a statistic that ranks how believable the 
event is to be an extragalactic MGF, with values for the true 
population is shown in orange. The background confidence 
intervals at 1, 3, and 50 are shown in blue. The four most 
significant events together surpass 5c discovery significance. 


Cline et al. 1982), which has even larger angular extent 
than M31, giving Q = 500. 

We perform a number of sanity checks to ensure our 
assumptions do not significantly affect our results. The 
search we run assuming our centroid a = 1.3 value; how- 
ever, we have confirmed that running the search at the 
90% confidence interval bounds (a = 0.5,2.2) identi- 
fies the same four bursts as significant outliers and does 
not identify other candidates. Running the search at 
greater NSIDE affects our Q values by <10%. Rerun- 
ning the search where the linear SFR weighting is al- 
tered to the stellar mass results in identification of the 
same galaxies but with generally lower Q values. Run- 
ning with specific SFR returns similar results. Together 
these suggest a progenitor that tracks SFR. Our results 
are insensitive to the assumed Fiso,min, so long as we do 
not exclude known events, as events of this strength are 
not detected far into the universe. There are a few events 
with €) > 1 which are either excluded as events of inter- 
est for our MGF search or insignificant given our sam- 
ple. Lastly, significantly raising the assumed Figo max 
marginally identifies GRB 100216A (Q = 10) which in- 
deed has a potential host galaxy within 200 Mpc (Perley 
et al. 2010), which is inconsistent with expectations for 
MGFs. 


3. PROGENITOR INVESTIGATIONS 


To determine the origin of these four bursts we first 
determine if the known GRB progenitors are compati- 
ble. Collapsars power long GRBs with durations 22s 
and are followed immediately by afterglow and then by 
broad-lined type Ic supernovae. This origin is excluded 
as all four events have durations 0.1s or less. Addition- 


ally, no subsequent supernova were reported in any case 
(Li et al. 2011b; though see Gehrels et al. 2006; Grupe 
et al. 2007). A neutron star merger origin is excluded by 
LIGO non-detections in gravitational waves for three of 
the four events (Abbott et al. 2008; Abadie et al. 2012; 
Aasi et al. 2014), but observations are insufficiently sen- 
sitive to inform on the origin of GRB 2004154. One may 
consider if off-axis GRBs could explain these events. 
The best known such event is GRB 170817A where the 
duration was longer and spectrum softer than the bulk 
of the short GRB population, which is inconsistent with 
the prompt emission from these four local events. Fur- 
ther, the rates of these local events (discussed in the 
following section) are orders of magnitude higher than 
cosmological GRBs (Siegel et al. 2019), even considering 
events that are oriented away from Earth. 

To determine the progenitors of these events we follow 
the historical procedure, where we begin by population 
comparison of prompt emission parameters. The only 
additional potential progenitor for extragalatic GRBs 
commonly discussed in the literature are MGFs where, 
contrary to the works that identified the two confirmed 
progenitors, we have the advantage of observations of 
galactic events which are summarized in Tablel. The 
parameters relevant for only the main peak of the flare 
that appear distinct from cosmological GRBs are the 
rise time and the intrinsic energetics. Figure4 contains 
the population comparison of these parameters. 

First, MGFs have rise times of order a few ms, far 
shorter than most cosmological short GRBs (Hakkila 
et al. 2018). Rise times are not reported in most GRB 
catalogs. As a proxy for the rise time we define the 
Time to Peak as the time from the start of the emission 
to the beginning of the peak 2 ms counts interval. An 
Anderson-Darling k-sample test against 75 bright Konus 
short GRBs (~15% brightest bursts detected by Konus 
between 1994 and 2020) rejects the null hypothesis that 
they are drawn from the same population at 299.996 
confidence. 

Second, MGF Zi,, values are orders of magnitude 
fainter than cosmological GRBs, where only the unusual 
GRB 170817A (Abbott et al. 2017) is comparable. This 
parameter depends on the distance to the source, which 
is not directly observable from prompt emission. For 
some cosmological GRBs direct distance (redshift) de- 
termination is made from follow-up observations. How- 
ever, for most short GRBs the distance is determined by 
first robustly associating the short GRB to an aligned or 
nearly aligned host galaxy, and then determining the dis- 
tance to the host (Fong et al. 2015). We adapt this last 
approach for MGFs to enable the use of larger prompt 
emission localizations and expected host galaxy proper- 
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Figure 3. The lightcurves of the candidate extragalactic MGFS in order of significance from Extended Data Table 1. These 
are from Konus- Wind and plotted with 2ms resolution (Frederiks et al. 2007b; Mazets et al. 2008; Svinkin et al. 2021), with 
GRB 070222 reported here for the first time. While GRBs 200415A and 051103 are strikingly similar (Svinkin et al. 2021) and 
GRB 070201 is broadly consistent with a single emission episode, GRB 070222 has two temporally and spectrally distinct pulses 


see Appendix B), suggesting varied behavior. 
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ties. For each GRB and potential host galaxy we cal- 
culate Qnost = 4r X; POP pHost/A; with PHost the 
weighted spatial distribution of that galaxy. Each GRB 
has only a single likely host, providing robust associa- 
tion. GRB 051103 has been discussed in the literature as 
belonging to the M81 Group of galaxies (Frederiks et al. 
2007b), which is dominating by the interacting galaxies 
M81 and M82. Our galaxy catalog selection and method 
assigns the burst to M82. 

The inferred Ej, values for each extragalatic MGF 
candidate is given in Table 1. For the population com- 
parison we add the Eiso distribution of GBM short 
GRBs (Abbott et al. 2017) to the sample of Konus bursts 
with measured redshift (Tsvetkova et al. 2017). To- 
gether these give 23 short GRBs with E;4 determined 
by a broadband instrument, which is the largest such 
sample to date. The extragalactic MGFs are clearly 
inconsistent with the broader population, rejecting the 
null hypothesis at 799.996 confidence. 


Host galaxy studies of GRBs have been key in de- 
termining prior progenitor channels (e.g. Fong et al. 
2015). As discussed in the design of our method, MGFs 
are expected to arise in star-forming galaxies or star- 
forming regions. Within our maximal detection distance 
for these bright events the galaxies with the highest 
SFR are M82, M83, NGC253, and NGC 4945 (Mattila 
et al. 2012). GRB 051103 is associated to M82 by our 
method or is consistent with star-forming knots on the 
outskirts of M81 (Ofek et al. 2006), GRB 070222 to M83, 
and GRB 2004154 to the star-forming core of NGC 253 
(Svinkin et al. 2021). GRB 790305B is associated to the 
star-forming Large Magellanic Cloud. This is consistent 
with a massive-star progenitor, as expected for an MGF 
origin. 

Individually, GRBs 200415A and 051103 are the most 
robust identifications of extragalactic MGFs based on 
our significance assessment and the results of part- 
ner analyses including lightcurve morphology and sub- 


Known Extragalactic 
MGF Event | 790305B 980827 041227 | 200415A 070222 051103 070201 
Origin 
False Alarm Rate 0 0 0 4.9x 107° 718x109 15x10? 1.2 x 1074 
BNS Excl. [Mpc] 6.7 5.2 3.5 
Galaxy Properties 
Catalog Name | LMC MW MW NGC253 M83 M82 M31 
Distance [Mpc 0.054 0.0125 0.0087 3.5 4.5 3.7 0.78 
SFR [Mo/yr 0.56 1.65 1.65 4.9 4.2 7.1 0.4 
GRB Properties 
Duration [s] | «0.25 «L0 «032 0.100 0.038 0.138 0.010 
Rise Time [ms ~2 ~4 ~1 2 4 2 24 
LYS? [10% erg/s 0.65 2.3 35 140 40 180 12 
Eiso [10% erg 0.7 0.43 23 13 6.2 53 1.6 
Index -0.7 0.0 -1.0 -0.2 -0.6 
Epeak [keV 500 1200 850 1080 1290 2150 280 


Table 1. A summary of the MGF sample. Significance for extragalactic events is from this text. BNS Excl. refers to the 
neutron star merger exclusion distances from LIGO. LMC refers to the Large Magellanic Cloud and MW refers to the Milky 
Way. Individual significance is determined by comparison of the individual Q against the full background sample. Distances 
for the known magnetars come from Olausen & Kaspi (2014); extragalactic distances are taken from the host galaxy values 
(which have minor variations with our catalog values). GRB parameters include Epeak as the energy of peak output, Index is 
the low-energy power-law from the spectral fit, and the rest are discussed in the text. GRB measures for the galactic events are 
from the literature; GRB measures for extragalactic events are all measured from Konus- wind data. 


millisecond variation of the prompt emission (Svinkin 
et al. 2021; Roberts et al. 2021). Newly identified is 
GRB070222 which is in-class with key properties of 
MGFs. However, it has two distinct but overlapping 
pulses, which is not known to occur from galactic events. 
This requires either a broader morphology of MGFs, a 
distinct and unknown origin, or a 1 in 100,000 chance 
alignment (Table 1). However, given the range of (quasi- 
)periodic oscillations seen from magnetar emission such 
a morphology is not necessarily surprising. 

To summarize the observational case for an MGF ori- 
gin: these events localize to the nearby universe and in 
particular to star-forming regions or star-forming galax- 
ies. The prompt emission is inconsistent with a collap- 
sar origin and gravitational wave observations exclude 
a compact merger involving neutron stars and/or black 
holes. The event rates, quantified below, are in excess 
of the majority of energetic astrophysical transients but 
are consistent with predictions from the known MGFs. 
The properties of the prompt emission are distinct from 
the larger short GRB population but again consistent 
with the properties from the known MGFs. There is 
additional evidence for individual events in partner anal- 
yses. We conclude that we have confirmed a sample of 


extragalatic MGFs that match prior predictions on de- 
tection rates and properties from both theoretical and 
observational studies. 

A remaining question is: why have we not identified 
MGFs to greater distances? Previously, MGFs were 
thought to be detectable to tens of Mpc. The spec- 
tra of the initial pulse of GRBs 200415A, 051103, and 
GRB 070222 are particularly spectrally hard with shal- 
low spectral index and high peak energies, which is con- 
sistent with GRB 041227 (Frederiks et al. 2007a). As- 
suming a cut-off power-law spectrum for bright MGFs 
with a low-energy spectral index = 0.0 and peak energies 
1.5 MeV produces only 15-20% of the photons in the 
nominal triggering energy range of 50-300 keV as com- 
pared to a typical short GRB (assuming an index of 0.4 
and peak energy 0.6 MeV; Goldstein et al. 2017). GRB 
monitors are triggered by photon counts, which sug- 
gests that the harder spectrum reduces the detectable 
by a factor of ^5 and therefore the volume by a fac- 
tor of more than 100. Instrument-specific comments are 
given in Appendix A. Further, there is a local overden- 
sity within +5 Mpc of CCSN (Mattila et al. 2012), which 
provides additional explanation of detections within this 
range and lack of detections beyond it. 
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Figure 4. Key parameter comparison of the extragalactic 
MGF candidates against the wider short GRB population 
and the known MGFs. Top shows the Time to Peak Konus 
distributions and bottom the Es distributions. The only 
comparable Fi value for a burst from a neutron star merger 
is the off-axis GRB 170817A. 


4. INFERENCES 


We now proceed to make population-level inferences 
utilizing the three known MGFs and treating all four of 
our events as extragalatic MGFs. 


4.1. Intrinsic Energetics Distribution 


The power-law distribution of the energetics of normal 
SGR flares gave hints to the physical process that pro- 
duces them (Cheng et al. 1996). Thus, it is interesting 
to measure the slope of the Fiso distribution for MGFs. 
We assign our search volume and detection threshold by 
empirical means, selecting 2.0 x 10-Sergem ? for the 
IPN and a maximal detection distance of +5 Mpc. We 
further restrict our sample to the past 27 years, where we 
have sufficient sensitivity to extragalactic events, leaving 
the 6 most recent bursts (excluding GRB 790305B). 

We assume the same power-law functional form for 
the Eiso PDF as our search method; however, we can- 


not utilize the maximum likelihood estimate because 
it requires the assumption that the observed sample is 
complete, which is not true for MGFs at extragalactic 
distances. Instead we simulate a large number of ex- 
tragalactic MGFs by drawing Eiso from PDFs over a 
range of a values, assigning them to specific host galax- 
ies weighted by their SFR, and setting the event dis- 
tance as the host galaxy distance. Events that would 
be detected are those where the sampled E;4 and dis- 
tance produce a flux greater than our detection thresh- 
old. Eisomin = 3.7 x 104+ erg is determined by sam- 
pling the Kolmogorov-Smirnov test statistic value over 
a range of viable options (Bauke 2007). Then, we cal- 
culate an Anderson-Darling k-sample value for a range 
of potentially viable a values. We take the 5% rejec- 
tion values as the bounds on a 9096 confidence interval, 
and determine the mean assuming a symmetric Gaus- 
sian distribution, giving a = 1.7 + 0.4. We note that 
this is consistent with the reported slope values of 5/3 
(Cheng et al. 1996) and 1.9 (Gótz et al. 2006) recurrent 
flares from galactic SGRs. 


4.2. Rates 


Utilizing the same sample and selection above we can 
constrain the intrinsic volumetric rate of MGFs. The 
dominant sources of uncertainty are the Poisson uncer- 
tainty and the imprecisely known sample completeness. 
The latter is limited by the uncertainty on the power- 
law index of the intrinsic energetics function, where for 
a steep index the majority of events will be missed 
(with most events below 1.0 x 10 erg missed in our 
sample volume) and for a shallow index most events 
are recovered. The o distribution is taken as a Gaus- 
sian. The SFR within 5 Mpc is 35.5 Mo/yr which is 
scaled to a volumetric rate by considering the total 
SFR within 50Mpc, which is 4000 M /yr from our 
galaxy sample. We infer a volumetric rate of Rygr = 
3.8*$9 x 105 Gpc 3 yr 1. 


4.3. Magnetar Formation Channel 


Magnetars may be generated in a variety of events 
including common CCSN, low-mass mergers (Price & 
Rosswog 2006), a rare evolution of white dwarfs (Dessart 
et al. 2007), or a rare sub-type of CCSN such as collap- 
sars or superluminous supernovae (Nicholl et al. 2017). 
Each of these is consistent with the observed association 
of magnetars to supernova remnants (Beniamini et al. 
2019). Low-mass merger events have long inspiral times 
and should track total stellar mass rather than the cur- 
rent SFR, which is disfavored given our model preference 
for SFR over stellar mass and the discovery of the first 
MGF from the LMC. A CCSN origin would arise from 
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regions with high rates of star formation. This is con- 
sistent with our observations and bolstered by both the 
lack of detections beyond 5 Mpc due to the local SFR 
overdensity and the detection of GRB 790305B from the 
low-mass, star-forming Large Magellanic Cloud. The 
host galaxies of our extragalactic sample and the Milky 
Way itself have larger mass and higher metallicity than 
is typically seen in hosts of collapsars or superluminous 
supernovae (Taggart & Perley 2019). "Therefore, the 
types of host galaxies favor common CCSN as the dom- 
inant formation channel of magnetars. 

Additional support for this conclusion is provided 
from the event rates. We can relate our inferred 
MGF rates to progenitor formation rates as RMGF = 
REventfMTActivermar/m (Tendulkar et al. 2016) where 
Revent is the rate of events that may form magnetars, 
fu is the fraction that successfully form magnetars, 
TActive the timescale that magnetars can produce MGFs, 
and rmar/m the rate of MGFs per magnetar. We take 
TActive © 104 yr limited by the decay of the magnetic 
field (Beniamini et al. 2019). Given the incompleteness 
of our known magnetar sample and lack of understand- 
ing which magnetars can produce MGFs, we use only 
the 3 known to be capable to estimate an upper bound 
of rugr;u < 0.02yr^! per magnetar. We note this is 
significantly weaker than those reported in the literature 
that consider all known SGRs, being ~ 1 x 1074 yr SGR 
(e.g. Ofek 2007; Svinkin et al. 2015). 

Of the discussed formation channels only CCSN are 
expected to track star-forming regions and have a com- 
parable rate, being 7x104 Gpc~* yr-! in the local uni- 
verse (Li et al. 2011a). A fiducial value on fm is 0.4 
with a 20 confidence interval of 0.12-1.0 (Beniamini 
et al. 2019); other estimates range between 0.01 and 
0.1 (e.g. Woods & Thompson 2004; Gullón et al. 2015). 
We require either that some magnetars produce multi- 
ple MGFs or that both fm œ% 1 and the true rate of 
Ruar is near our 95% lower bound. Alternatively, us- 
ing the CCSN rate and the 95% lower limit on Rar 
we can place observational constraints using our results 
of fm > 0.005, further excluding particularly rare sub- 
types of, and favoring common, CCSN as the dominant 
formation channel of magnetars. 


5. CONCLUSIONS 
To summarize our conclusions: 
e We have shown that 4 short GRBs occurred within 
~5 Mpc which are the closest events by an order of 


magnitude in distance. Our analysis was the first 
to identify GRB 070222 as a local event. 


e They are inconsistent with a collapsar or neutron 
star merger origin. 


e Their prompt emission is inconsistent with the 
properties of cosmological GRBs, but is consistent 
with the observations of the known MGFs. 


They originate from star-forming regions or star- 
forming galaxies, including those with metallicity 
that prevents collapsars from occurring. 


Altogether this matches expectations for an MGF 
origin, which appear to produce 4 out of 250 
events. This would be ~2% of detected short 
GRBs (consistent with the 1-896 range from the 
literature Ofek 2007; Svinkin et al. 2015) or ~0.3% 
of all detected GRBs. 


Modeling the intrinsic energetics distribution of 
MGFs as a power-law constrains the index to be 
1.7+ 0.4. 


The volumetric rates are Rygr = Bt x 
10° Gpe~3 yr-!. 


The rates and host galaxies of these events favor 
CCSN as the dominant formation channel for mag- 
netars, requiring at least 0.5% of CCSN to produce 
magnetars. 


e We estimate the rate of MGFs per magnetar to be 
< 0.02 yr-!. 


Our results suggest that some magnetars produce 
multiple MGFs: this would be the first known 
source of repeating GRBs. 


GRB 070222 suggests MGFs can have multiple 
pulses. 


e MGFs may not be detectable to tens of Mpc with 
existing instruments due to their spectral hard- 
ness. 


Our analysis suggests additional extragalactic MGFs 
may be identified with improved analysis but *smoking- 
gun" confirmation likely requires future instruments. 
The inferred rates are sufficiently high that they may 
contribute to the stochastic background of gravitational 
waves. This, and the recent observations of a fast radio 
burst to lower-energy gamma-ray flares from magnetars 
(Bochenek et al. 2020; Marcote et al. 2020; Ridnaia et al. 
2020; Li et al. 2020), suggest the coming years will bring 
new insights into the physics and emission of magnetars. 
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6. APPENDIX A 


We present rough estimates for the maximal detection distance of bright MGFs with representative active instru- 
ments. Konus- Wind can detect bright MGFs to ~13-16 Mpc, based on GRBs 051103 and 200415A (Svinkin et al. 
2021). This can be taken as the approximate detection distance of the IPN (Svinkin et al. 2015). The following inves- 
tigations assume a hard spectrum based on the time-integrated values for the most energetic bursts, with a low-energy 
spectral index ~ 0.0 and peak energies &1.5 MeV. This has only 15% (20%) the number of photons over the 15-150 keV 
(50-300 keV) energy range, reducing the detection distance by ~x5 and thus a reduction in volume of 7100. 

The GBM GRB trigger algorithms cover 50-300 keV where the short GRB sensitivity is usually quoted over the 
64ms timescale. With the assumed spectral and energetics values GBM would have only triggered these on-board 
algorithms out to ~15-20 Mpc. At greater distances only the peak flux interval would be visible, which would be 
spectrally harder, and reduce this distance. GBM localizations alone are insufficient to associate events to any specific 
burst. Ground-based searches for GRBs and Terrestrial Gamma-ray Flashes should be able to recover additional 
events, but may require confirmation in other GRB instruments 

The INTEGRAL SPI-ACS and IBIS are especially sensitive to hard and short bursts, and additional extragalactic 
MGFs have likely triggered the SPI-ACS real-time pipeline in the past. However, SPI-ACS and IBIS lack the capacity to 
discriminate extragalactic MGFs from high cosmic ray effects which appear similar to real events in these instruments. 
The real-time IBAS pipeline has not been tuned to favor short and hard events. We estimate that SPI-ACS would 
record sufficient signal from extragalactic MGFs for association with another instrument up to 25-35 Mpc, but would 
only independently report much brighter events out to 15-20 Mpc. Sensitivity of IBIS is close or better than to that 
of SPI-ACS in about 1096 of the sky, and in the majority of directions, IBIS would only yield detectable signal for 
extragalactic MGF flares out to at most 10 Mpc. However, PICsIT may often be more suitable for triangulation, 
owing to better time resolution, and can provide some spectral characterization. 

The Swift BAT has >500 different rate trigger criteria running in real-time onboard, continuously sampling and 
testing trigger timescales from 4ms up to 64 seconds, each of which is evaluated for 36 different combinations of energy 
ranges and focal plane regions. While the BAT detector is sensitive to photons with energies up to 500 keV, the 
transparency of the lead tiles in the mask above 200 keV limits its imaging energy range (necessary for a successful 
autonomous trigger) to 15-150 keV. This narrow and low energy range limits the BAT’s sensitivity to hard events, 
such as MGFs, despite its high effective area. Due to the number and complexity of the onboard triggering algorithms, 
the varying compute load on the BAT CPU, as well as the evolving state of the BAT detector array and changing 
operational choices for trigger vetoes/thresholds, modelling the likelihood of an onboard autonomous trigger is quite 
difficult. In addition, due to BAT's high effective area, continuous time-tagged event data cannot be downlinked, 
making it difficult to assess the relative completeness of the triggering algorithms vs ground searches, though this 
is partly ameliorated by GUANO (Tohuvavohu et al. 2020). Under the assumed energetics and spectral values, we 
estimate that as of 2020 (averaging half of the original detector array online) Swift/BAT should reliably trigger on 
MGFs out to ~25 Mpc in the highest coded region of its field of view. Ground analyses in the downlinked BAT 
event data can extend this, but the availability of this data will often depend on an external trigger (e.g. GUANO). 
We note that operational changes to the BAT onboard triggering thresholds with the goal of increasing sensitivity 
to extragalactic MGFs and local low-luminosity GRBs have been previously attempted. In 2012 the threshold for a 
successful trigger from an image was lowered from the usual value of 6.5 to 5.7, with the condition that triggers in this 
range be localized to within 12 arcminute projected offset from a local catalogued galaxy stored in the BAT onboard 
catalog. No local GRB-like source was ever identified in this program. 


7. APPENDIX B 


As GRB 070222 has not been reported elsewhere we describe its basic analysis here. The event was detected by 
Konus- Wind, HEND on Mars Odyssey, and both SPI-ACS and PICsIT on INTEGRAL. Combination of the two best 
annuli produce a localization with a 90% containment region of 0.004 deg?. This location and its consistency with M83 
is shown in Figure 5. 

'This burst is distinct from the separate candidates as having two separate pulses. Time-resolved analysis of this 
burst is summarized in Table2 while time integrated analysis is reported in the Second Konus GRB Catalog (Svinkin 
et al. 2016). 
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Figure 5. The localization of GRB 070222 compared to the position and angular size of M83. 


Tstart | Tstop Index Epeak Flux 

[s] [s] [keV] | [1 x 1078 erg/s/cm?] 
-0.006 | 0.012 | 0.144923 | 7331335 153.4* 258 
0.026 | 0.038 | —0.27*035 | 193135 24.573 5 


Table 2. The time-resolved analysis of the two pulses of GRB 070222. Errors are quoted at 68% confidence. The main pulse 
is spectrally hard, similar to the time-integrated fits of GRB 200415A and GRB 051103. 


